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ABSTRACT 

The thermal history of the intergalactic medium (IGM) after reionization is to a large extent deter- 
mined by photoheating. Here we demonstrate that calculations of the photoheating rate which neglect 
radiative transfer effects substantially underestimate the energy input during and after reionization. The 
neglect of radiative transfer effects results in temperatures of the IGM which are too low by a factor 
of two after He n reionization. We briefly discuss implications for the absorption properties of the IGM 
and the distribution of baryons in shallow potential wells. 

Subject headings: quasars: absorption lines - galaxies: formation - cosmology: theory 
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1. INTRODUCTION 

The absence of a Gunn- Peterson trough (Gunn & Peter- 
son 1965, Scheuer 1965) in the spectra of high-redshift ob- 
jects is solid evidence that the Universe has been reionized 
before z ~ 5. Hydrodynamical simulations have shown 
convincingly that at redshifts z ~ 3 most baryons are still 
contained in a photoionized IGM which is responsible for 
the Lya forest in the absorption spectra of high-redshift 
QSO's (Cen et al. 1994, Petitjean, Mucket & Kates 1995; 
Zhang, Anninos & Norman 1995; Hernquist et al. 1996). 
However, the detailed distribution of the baryonic compo- 
nent as well as the absorption properties of the IGM will 
depend on the thermal history of the IGM and thus the 
heat input due to photoheating (Efstathiou 1992, Miralda- 
Escude & Rees 1994, Navarro & Steinmetz 1997, Haehnelt 
& Steinmetz 1998, Bryan et al. 1999; Theuns et al. 1999). 
Most numerical simulations so far adopted the optically 
thin limit when calculating photoheating rates. This is, 
however, a bad approximation during reionization when 
the optical depth is large. As we will argue later in 
the case of helium it is also a bad approximation after 
He ii reionization (see also Giroux & Shapiro 1996). Here 
we investigate radiative transfer effects during and after 
reionization and discuss implications for the thermal his- 
tory and absorption properties of the IGM and the distri- 
bution of baryons in shallow potential wells. We assume 
an Einstein-de-Sitter Universe and a Hubble constant of 
50kms _1 Mpc" 1 . 



2.1. 



2. PHOTOHEATING 
Optical thin vs optical thick 



In the optically thin limit the mean excess energy of 
ionizing photons is given by 
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where J(y) is the flux of ionizing photons o(y) is the 
crossection for ionization, v t h the threshold energy for 
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ionization, and h Planck's constant. For a powerlaw 
spectrum J{v) = Jo (v/vth)~ a and a powerlaw approxi- 
mation (oc v~~ 3 ) for the crossection one finds (E^ n ) w 
ft,z/ t h/(a + 2) = E t h/(a + 2). For spectral indices of in- 
terest a > 1 these excess energies are small fractions of 
the threshold energy E t h- This is because the integral 
in equation (1) is strongly weighted towards the threshold 
energy. In the optical thick limit where every ionizing pho- 
ton emitted is absorbed the mean excess energy is typically 
much larger than (i?*^ 111 ) of equation (1). For the power- 
law spectrum discussed above (E^ ck ) w hv t h/(ct — 1). 

2.2. Energy input during reionization 

A non-equilibrium multi-frequency 1-D radiative trans- 
fer code is used to study the thermal structure of R-type 
ionization fronts surrounding a point source in spherical 
symmetry (Abel, Meiksin, & Norman 1999). The code 
evolves the cooling/heating and chemistry model for pri- 
mordial gas of Abel et al. (1997) with the methods dis- 
cussed by Anninos et al. (1997). It solves the static radia- 
tive transfer equation in an expanding universe and explic- 
itly conserves photons (Abel, Norman and Madau 1999). 

Figure 1 shows radial temperature profiles of two dif- 
ferent ionization fronts around a luminous QSO (photon 
emission rate A ph = 3 x 10 56 s" 1 , e(v) = ea [v / vq)~ xa ) 
at z = 6 and z = 4, respectively. In both cases the front 
evolves into a homogeneous IGM of primordial composi- 
tion and mean baryonic density (fl^h 2 — 0.02). In the 
first case the IGM is assumed to be completely neutral 
and hydrogen and helium reionization occur almost simul- 
taneously. In the second case the gas is assumed to be 
predominantly in the form of HII and Hell due to a soft 
UV background already present when the QSO switches 
on. For both models the ionization front is shown after 
10 7 and 10 8 yr. 

In the first model the ionization front is very sharp due 
to the small mean free path of hydrogen ionizing pho- 
tons and the temperature behind the reionization front 
rises to 50000 K. This corresponds to (i? p h) ~ 0.85 x _E t h- 
This is much larger than in the optical thin approximation, 
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(^ph' n ) = 0.26^th, and approaches the optical thick limit, 
(^ph lck ) = 1-25-Eth, for the adopted power-law spectrum 
with a = 1.8. The difference is mainly due to line-cooling 
of the neutral component in the ionization front. This be- 
comes more apparent at later times when the ionization 
front has slowed down. Hydrogen line cooling of the neu- 
tral component during the passage of the ionization front 
becomes then more effective and the temperature falls to 
28000 K (Miralda-Escude & Rees 1994). This is still nearly 
a factor two larger than the value of 15000i4T in the optical 
thin limit (e.g. Haehnelt & Steinmetz 1998). 
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Fig. 1. — Thin and thicS curves show Ino radial temperature 
profiles of the ionized region around a luminous quasar (N p ^ = 
3 X 10 56 s _1 ) that switched on at z = 6 and z = 4, respectively. 
In the first case the ionization front evolves into a completely neu- 
tral homogeneous IGM of primordial composition while in the latter 
case the gas is assumed to be predominatly in the form of H II and 
He II . For the second case the inset shows the spectral shape of the 
ionizing flux at the three different radii marked in the figure. 

In the second case only helium is reionized and hydrogen 
line cooling plays no role due to the small neutral fraction 
of hydrogen. The mean free path of Hell ionizing photons 
is about a factor fifty larger than that of hydrogen ionizing 
photons and the ionization front extends now over a few 
(comoving) Mpc. The temperature is raised from 9000K 
to about 17000K and the total energy input corresponds to 
(E p h) ~ 0.5 x Eth, but this value rises to {E p h) ~ 0.7 x E t h 
when the energy deposited ahead of the front is taken into 
account. The value is again smaller than that in the op- 
tical thick limit. This time because the mean free path of 
an helium ionizing photon becomes larger than the size of 
the ionized region for (1 + E ph /E th ) > 3.8(i?/30 Mpc) 1 / 3 , 
where we have assumed that the ionization cross section 
scales as v~ z above the threshold. This also explains why 
the IGM temperature behind the ionization front increases 
with increasing radius. The inset in Figure 1 shows the 
spectral shape behind, within and ahead of the front to 
demonstrate how the spectrum of absorbed photons be- 
comes harder across the ionization front. 

2.3. Photoheating after reonization 

After reionization the mean free path of hydrogen ion- 
izing photons is comparable to the mean distance be- 
tween hydrogen Lyman limit systems (LLS) which is about 
350((1 + z)/4)~ 3 Mpc comoving (Storrie-Lombardi 1994). 



Much less is known about the number of helium LLS. 
The model of Miralda-Escude, Haehnelt, & Rees (1999, 
MHR99) which uses a realistic density distribution to treat 
inhomogeneous reionization predicts a typical distance of 
~ 30 Mpc around z ~ 3 (see also Fardal, Giroux & Shull 
1998). Both values are larger than the typical size of nu- 
merical simulations and this is is often taken as argument 
to treat photoheating in the optically thin limit. However, 
it is not before z=1.6 that the mean free path of hydrogen 
ionizing photons becomes comparable to the Hubble ra- 
dius (e.g. Madau, Haardt & Rees 1999) and the Universe 
is likely to be opaque to helium ionizing photons up to the 
present day. Almost all photons emitted at high redshift 
will therefore be absorbed within a Hubble time and cal- 
culations adopting the optically thin limit must severely 
underestimate at least the mass-weighted average photo- 
heating rate after reionization. 

In the case of hydrogen the mean free path of an hydro- 
gen ionizing photon is smaller than the distance between 
ionizing sources and the optical depth of LLS is large. At 
a typical place in the Universe the UV background will be 
dominated by ionizing photons emitted within the region 
filling the volume to the next LLSs. In the case of Hell 
the mean distance between sources of ionizing photons at 
z ~ 3 is ~ 100 Mpc (MHR99) and thus still larger than 
the mean distance between helium LLS. At a typical point 
in the Universe the Hell ionizing photons will therefore- 
have passed a helium LLS. The minimum energy input per 
ionization is thus raised such that (1 + E p h/E t h) > t 1 / 3 , 
where r is the optical depth at the ionization threshold. 
Little is known about the optical depth distribution of He- 
lium LLS but even a rather low optical depth of r ~ 10 
would raise the Hell photoheating rate by more than a 
factor of four compared to the optical thin case. 

However, even when the distance between helium LLS 
falls below the mean distance between sources of ionizing 
photons the optical thin approximation will significantly 
underestimate the Hell photo-heating rate, unless helium 
LLS have a very large optical depth. To understand this 
consider the extreme case of a homogeneous emissivity of 
ionizing photons and many helium LLS with optical depth 
of order unity at the He n ionization threshold. The ioniz- 
ing background at a given frequency would then be given 
by the emisssivity times the mean free path at this fre- 
quency, I{v) oc e{v)l(v) cx e(v)<j(v)~ l . In this limit the 
mean absorbed photon energy is equal to the mean emit- 
ted energy of photons. 

An accurate prediction of the mean absorbed energy is 
extremely difficult as it depends very sensitively on the 
spatial distribution of Hell optical depths and ionizing 
sources during and after reionization and on their emit- 
ted spectrum. We can here only roughly estimate that 
Eph/ Eth ~ 0.5 to 1 at z ~ 3 and should decrease with 
redshift. The effective Hell photoheating rate should thus 
be between 2 and 4 times larger than that in the optical 
thin approximation. It could be even larger, if the emitted 
spectrum were harder than the a = 1.8 power law assumed 
here. 
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Fig. 2. — a) The temperature evolution of an IGM of primordial 
composition and mean density assuming a homogeneous UV back- 
ground as in equation (3). For the solid curves the Hell photoheat- 
ing rate was raised by a factor 2,3,4, respectively. The dashed curves 
assume a homogeneous UV background as given by Haardt & Madau 
(1996) with the Hell photohcating for the optical thin limit and 
raised by a factor four, respectively, b) The slope /3 = dlnT/dlnp 
of the temperature density relation at mean density. 

3. IMPLICATIONS 

3.1. Temperature of the IGM 

In the last section we have demonstrated that optical 
depth effects raise the effective helium photoheating rate 
by about a factor 2 to 4. Here we investigate the effect of 
such an increased helium photoheating rate on the tem- 
perature of the IGM. We used a non-equilibrium code to 
calculate the redshift evolution of the ionization state and 
temperature of an photoionized IGM of primordial com- 
position with mean baryonic density. We have assumed a 
homogeneous UV background of the form, 
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This is similar to model C in Haehnelt & Steinmetz 1998. 
The first term mimicks reionization of hydrogen by stars 
before z = 6. The second component represents the contri- 
bution of QSO's to the UV background which here results 
in a separate epoch of Hell reionization around z = 4. 

The dotted curves shows the temperature evolution for 
the optical thin case. For the solid curves we have raised 
the Hell photoheating rate by factors 2,3 and 4. This 
corresponds to an increase in the total photoheating rate 
by a factor 1.3, 2 and 2.6, respectively. For reference the 
dashed curves show the temperature evolution for the UV 
background as calculated by Haardt & Madau (1996) in 
the optical thin limit and with the Hell photohcating rate 
increased by a factor of four. 

As pointed out by by MHR99 the assumption of an 
homogeneous UV background during reionization results 
in an artificially abrupt reionization. In reality reioniza- 
tion should occur over a significant fraction of the Hub- 



ble time. The assumption of a homogeneous background 
during reionization should nevertheless give a resonable 
estimate of the mean temperature after reionization even 
though there will be considerable spatial temperature fluc- 
tuations. 

At z < 6 adiabatic cooling exceeds Compton cooling and 
after helium reionization and in a fully photoionized gas 
the temperature is thus determined by the balance be- 
tween the photoheating rate, r p h = (aHeiii(-E'p i h III )^Heiii + 
aHii(-E ; ™ I )wHii)^e and the adiabatic cooling rate C a( j = 
3kTH(z)p/ firrip, where a(oc T~ ' 7 ) is the recombination 
coefficient , n e is the electron density, H{z) is the Hubble 
constant /i is the mean molecular weight and m p is the 
proton mass. Note that there is no explicit dependence on 
the strength of the radiation field. This results in a quasi- 
equilibrium temperature which is approached at late times 
(Miralda Escude & Rees 1994) 
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where we have taken the optical thin limit for hydrogen 
heating. The increased Hell photoheating rate due to op- 
tical depth effects results in an increase of the temperature 
by factors 1.5-2.5 (Fig. 2a). Changing to a cosmological 
model with £1 = 0.3 would further raise this temperature 
by 20 percent at redshift three because of the decrease in 
the adiabatic cooling rate. 

Madau & Efsthathiou (1999) recently suggested that 
Compton heating by the hard X-ray background might 
significantly contribute to the heating of the IGM. They 
obtain for the Compton heating rate at the present day 
H co = 1.25 x 10~ 31 erg s _1 . We can now compare that to 
the photoheating rate 
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Note that the photoheating rate incrases with density 
while the Compton heating rate does not. The redshift 
evolution of the Compton heating rate f(z) depends on 
the evolution of the energy density of the hard X-ray back- 
ground the mean energy of the Compton scattered photon. 
The Klein-Nishina reduction to the Compton cross section 
becomes increasingly important at high redshift. Madau & 
Eftathiou (1999) obtain f(z) ~ (l + z) 13 / 3 as approximate 
redshift dependence if the hard X-ray background was pro- 
duced at high redshift. Hence, the Compton heating rate 
by the hard X-ray background exceeds the photo-heating 
rate only at redshifts larger than 
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3.2. Absorption properties of the IGM 

The raise in temperature should affect the absorption 
properties of the IGM in several ways. It will increase 
the thermal broadening of absorption lines but might 
also introduce additional non-thermal broadening due to 
pressure-induced motions. The temperature increase due 
to optical depth effects during Hell photohcating is there- 
fore a good candidate to resolve recently reported diffi- 
culties of numerical simulations to reproduce Dopplcr pa- 
rameters of Lya absorption lines at z — 2 to 3 as large as 
those observed (Bryan et al. 1999; Theuns et al. 1999). 
This would otherwise require an unreasonable large value 
of fl^h 2 (Theuns et al 1999). The temperature increase 
will also increase the mean baryonic density required to 
match the mean flux decrement in QSO absorption spec- 
tra (Rauch et al. 1997). The details of helium heating will 
furthermore affect the temperature density relation (Hui 
& Gnedin 1997) and thus the flux decrement distribution. 
To demonstrate the effect on the temperature density re- 
lation we have plotted in Fig. 2b the logarithmic slope 
ft = d In T I d In p at mean density for the same models as 
in Fig 2a. 

3.3. Baryons in shallow potential wells 

The increased helium photoheating rate will also raise 
the Jeans mass and will decrease the ability of the gas to 
cool and collapse into shallow potential wells. It will thus 
affect the ill distribution of baryons. For photoion- 
ized gas the dominant cooling processes during the early 
stages of the collapse are recombination cooling and Brem- 
strahlung. The equilibrium between photohcating and re- 
combination cooling (C roc oc p 2 T 3 ) results in a temper- 
ature T roc ~ 2.3 x 10 5 [{(E^ ul )/E^ cm ) + 0.15]K. The 
increased Hell photoheating rate should thus affect the cir- 
cular velocity threshold above which gas can collapse into 
dark matter haloes and the fraction of gas that can cool in 
shallow potential wells (Thoul& Weinberg 1996, Forcada- 
Miro 1998). Similarly the increased hydrogen and he- 
lium photoheating rates should raise the circular velocity 
threshold below which baryons in shallow potential wells 
which collapsed before reionization can be evaporated by 
an ionization front (Shapiro et al. 1998; Abel & Mo 1998; 
Barkana & Loeb 1999). 



4. CONCLUSIONS 

We have studied radiative transfer effects on the energy 
input into the IGM due to the reionization and photoheat- 
ing of hydrogen and helium. The energy input depends 
crucially on the spectral distribution of the photons actu- 
ally absorbed. 

During reionization the spectrum of an ionizing source 
becomes harder across the ionization front. The energy 
input is close to the mean energy of the photons between 
the ionization threshold and the frequency for which the 
mean free path of an ionizing photons becomes larger than 
the radius of the ionized region. 

After the reionization of Hell the mean spectrum of ion- 
izing photons above the He in ionization edge will also be 
considerably harder than the emitted spectrum of an in- 
dividual source. This is because hard photons can pass 
through helium LLS. The number and optical depth dis- 
tribution of helium LLS is still highly uncertain. Thus 
their effect on the UV specrtum is likely to have been un- 
derestimated previously. 

The increased Hell photohcating rate due to optical 
depth effects raises the temperature of the IGM typically 
by about a factor 1.5 to 2.5 compared to the optical thin 
approximation. It might thus resolve the problem that cur- 
rent numerical simulations of the Lya forest (which treat 
photoheating in the optical thin limit) produce Doppler 
parameters narrower than observed. 

The increased Hell photoheating rate will reduce the 
ability of the gas to collapse into shallow potential wells. 
It will also raise the circular velocity threshold below which 
gas in dark matter haloes is photoevaporated during reion- 
ization. 
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